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ABSTRACT 

Upcoming weak lensing surveys on wide fields will provide the opportunity to reconstruct the structure 
along the line of sight tomographically by employing photometric redshift information about the source 
distribution. We define power-spectrum statistics, including cross correlation between redshift bins, 
quantify the improvement that redshift information can make in cosmological parameter estimation, and 
discuss ways to optimize the redshift binning. We find that within the adiabatic cold dark matter class of 
models, crude tomography using two or three redshift bins is sufficient to extract most of the information 
and improve the measurements of cosmological parameters that determine the growth rate of structure 
by up to an order of magnitude. 

Subject headings: cosmology: theory - gravitational lensing - large-scale structure of universe 



1. INTRODUCTION 
With new instruments such as MEGACAM at CFHT 



( |Boulade et aL 19 98) and the VST a t the European South- 



ern Observatory (Arnaboldi et al. 1998), wide-field sur 



veys detect the weak lensing of faint galax ies by large 
scale structure will soon become a reality (see Melier 1998 
for a recent review). Weak lensing by large-scale struc- 
ture produces a correlated distorti on in the clliptici ties of 



the galaxies on the p ercent l evel (Blandford et al. 1991; 



Miralda-Escude 1991; Kaiser 1992) which can be used to 



measure a two-di mensional projection of the in tervening 
mass di stribution (Tyson, Valdes & Wenk 1990; Kaiser & 



Squires 1993) 



If the redshift of the source galaxies are known, then 
more information can be extracted out of weak lensing 
by tomography, i.e. differencing the two-dimensional pro- 
jected images to recover the three-dimensional distribu- 
tion. In the absence of spectroscopy, approximate red- 
shifts for the faint galaxies ca n be determ ined through 



photometric techniques (see e.g. Hogg et al. 1998 and ref- 



erences therein) and with the large number of galaxies at 
R <J 25 (~ 10 5 deg -2 ) the properties of distribution can be 
known to good accuracy (Seljak 1998). Indeed the weak 
lensing surveys already plan to use photometric redshift 
information at least on a small subsample to measure the 
low order moments of the distribution such as its mean. 
These are important f or determinin g the cosmologic al im- 



plications of the d ata (Smail et al. 1995; Fort et al. 1995; 



Luppino & Kaiser 1997) 



The potential of tomographic techniques, especially in 
the wide-field limit where the cosmological information is 
completely contained in the two-point functions or power 
spectra, remains larg ely unexplored . Ind eed mo st stu dies 
of weak lens ing (e.g. Jain & Seljak 1997; Kaiser 1998; |Ht 



k Tcgmark 1999) simply" assume a delta-function distri- 
bution of galaxies making tomography impossible. 

In this Letter, we study the prospects for weak lensing 
tomography within the framework of the adiabatic cold 
dark matter (CDM) class of models for structure forma- 
tion. We begin by defining the power spectrum statistics 
for an arbitrary set of galaxy redshift distributions. These 
are the power spectrum of the convergence map for each 



distribution and the cross-correlation between the maps. 
We then quantify how much additional information can be 
extracted by subdividing a single magnitude-limited sam- 
ple into bins in redshift and analyzing their joint power 
spectra and cross-correlation. We conclude with a discus- 
sion of how errors in photometric redshifts might affect 
tomographic techniques. 

2. POWER SPECTRA 



Generalizing the results of Kaiser (1992, 1998), we de- 
fine the angular power spectra and cross-correlation of sky 
maps of the convergence based o n a ser ies of galaxy red- 
shift distributions rii(z) (see also Seljak 1998) 
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where a™ m -j j are the spherical harmonic coefficients of the 
maps. HereD = J*(H /H)d z is the dimensionless comov- 
ing distance and 



d a {d) = n~ 



(2) 



is the dimensionless angular diameter distance, where 
k = 1 — J2i i s ^e effective density in spatial cur- 
vature in units of the critical density. The efficiency 
with which gravitational potential fluctuations as mea- 
sured by their dimensionless power per logarithmic interval 
A| = k 3 P$/2n 2 , lens the given galaxy distribution rii is 
described by 
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Note that rij is normalized so that J °° dzrii(z) = 1. Fi- 
nally kg = £Ho/ Da(D) is the wavenumber which projects 
onto the angular scale I at distance D. For small fields of 
view, the spherical harmonics of order I can be replaced 
by Fourier m odes with angular frequency w. 



We use the 



Peacock fc Dodds (1996) scaling relations to 



evalute A| in the non-linear density regime. Equation ([!]) 
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Fig. 1. — Subdividing the source population. Partitioning the 
galaxies by the median redshift (or distance D) yields lensing effi- 
ciencies with strong overlap. 



assumes that the redshift distributions are sufficiently wide 
to encompass many wavelengths of the relevant fluctua- 
tions (2w/ke) along the line of sight so that th e Limber 
equation holds even tomographically (see Kaiserj 1998). 

These power spectra define the cosmic signal. Shot noise 
in the measurement from the intrinsic ellipticity of the 
galaxies adds white noise to the cosmic signal making the 
observed power spectra 



c ij (i) = ^ j (£) + ( 1 l t )s ij /n i 



(4) 



where ( 7 ? t > J 
nent, 



yint) 1 2 i s the rms intrinsic shear in each compo- 
and fii is the number density of the galaxies per 
steradian on the sky in the whole distribution rii{z). 

The distributions rii(z) need not be physically distinct 
galaxy populations. Consider a total distribution n(z) 
with 

dz 1 

i— {D)<xD a vq>[-{D/D*Y\, (5) 



which roughly approximates that of a magnitude-limited 
survey, and take a = 1,(3 = 4 for definiteness (assumed 
throughout unless otherwise stated). One can subdivide 
the sample into redshift bins to define the distributions 
rii(z). The power spectra for cruder partitions can always 
be constructed out of finer ones: if the j and k bins are 
combined, then 
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In Fig. [T| we show an example where the galaxies with 
z < Zmedian are binned into n\ and the rest into ri2- Here 
and throughout we will take our fiducial cosmology as an 
adiabatic CDM model with matter density il m = 0.35, 
dimensionless Hubble constant h = 0.65, baryon density 
SI5 = 0.05, cosmological constant f^A = 0.65, neutrino 
mass m„ — 0.7 eV, the initial potential power spectrum 
amplitude A, and tilt ns = 1. 

We also plot in Fig. |l| the lensing efficiency func- 
tion gi{D). Notice that despite having non-overlapping 




Fig. 2. — Power spectra and cross correlation for a subdivision in 
two across the median redshift ^ me dian = 1 an d errors for a survey 

of 5° on the side, (7 i 2 nt ) 1/2 = 0.4, and n = 2 X 10 5 deg~ 2 . Note 
the strong correlation R^j between the two power spectra make the 
combination of the power spectra less constraining than a naive 
interpretation of the individual errors would imply. 



source distributions (upper panel), the lensing efficiencies 
strongly overlap (bottom panel) implying that the result- 
ing convergence maps will have a correspondingly large 
cross correlation. This is of course because the high and 
low redshift galaxies alike are lensed by low-redshift struc- 
tures. Also for this reason, there will be always be a 
stronger signal in the high redshift bins. This fact will 
be important for signal-to-noise considerations in choos- 
ing the bins. 

All of these properties can be seen in Fig. |^, where we 
plot the resultant power spectra and their cross correlation 
for the equal binning of Fig. pi 

3. REDSHIFT BINNING AND PARAMETER ESTIMATION 

While subdividing the sample into finer bins always in- 
creases the amount of information, there are two consid- 
erations that limit the effectiveness of redshift divisions. 
The first is set by the shot noise from the intrinsic ellip- 
ticities of the galaxies. Once the number density fii per 
bin is so small that shot noise surpasses the signal in equa- 
tion (Q) , further subdivision no longer helps. The point at 
which this occurs depends on the angular scale of inter- 
est. The greater number of galaxies encompassed by the 
larg er angu lar scales boosts the signal to noise (see Fig. || 
and Kaisei 1992). Based on this criterion, one should sep- 
arately subdivide the data to extract the maximal large 
and small angle information. 

However there is a second consideration. If the lens- 
ing signal does not change significantly across the red- 
shift range of the whole distribution, then subdivision will 
not add information. These considerations can be quanti- 
fied by considering the correlation coefficient between the 
power spectra of the subdivisions: Ry = P^ ) /(P^P" J ) 1 / 2 . 
For the model of Fig. 0, the power spectra are highly cor- 
related (R12 ~ 0.8) even with only two subdivisions. Thus 
even though there is enough signal to noise to subdivide 
the sample further, one gains little information by doing 
so. 
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One can combine these two considerations by diagonal- 
izing the covariance matrix and considering the signal to 
noise in the diagonal basis. The appropriate strategy for 
subdivision depends on the true redshift distribution of 
the galaxies and the model for structure formation. One 
should therefore perform this test on the actual data to 
decide how to subdivide the sample. 

Nevertheless, to make these considerations more con- 
crete, let us consider the specific goal of measuring the 
cosmological parameters p a assuming that the underlying 
adiabatic CDM cosmology described above is correct. The 
Fisher information matrix can be used to quantify the ef- 
fect of subdivision. It is defined as 
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where L is the likelihood of observing a data set x given 

the true parameters pi . . . p a - 

Generalizing the results of Hu fc Tegmark (1998) to mul- 
tiple correlated power spectra, we obtain[J] 




0.4 0.6 
upper fraction 

Fig. 3. — Tomographic error improvements on Cl\ for ^median = 1- 
Upper panel: improvement as a function of the fraction of galaxies 
in the upper redshift bin for 2 bins versus 3 bins (same fraction in 
upper two bins). Lower panel: redshift corresponding to the upper 
division. 
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under the assumption of Gaussian signal and noise, where 
/sky is fraction of sky covered by the survey, the covari- 
ance matrix C was defined in equation (|j), and commas 
denote partial derivatives with respect to the cosmological 
parameters p a . We take ^ max = 3000 to approximate the 
increased covariance due to th e nonlinearities producing 
non-Gaussianity in the signal ( [Scoccimarro, Zaldarriaga 
& Hu| 1999). Since the variance of an unbiased estimator 
of a parameter p a cannot be less than a(p a ) = (F _1 ) QQ , 
the Fisher matrix quantifies the best statistical errors on 
parameters possible with a given data set. 

For the purposes of this work, the absolute errors on 
parameters are less relevant than the improvement in er- 
rors from subdividing the data (see Hu fc Tegmark for 
the former). We therefore test a 4 dimensional subset of 
the adiabatic CDM parameter space to see how subdivi- 
sion can help separate initial power (A) from the various 
contributors to the redshift-dependent evolution of power 
(Q\,Qk ,m u ). For reference the standard errors a a for this 
parameter space without subdivision are given in Table 
1. Errors in the full parameter space would be increased 
but note that the neglected parameters (£l m h 2 ,tti,h 2 , and 
ns) are exactly those that the CM B satellite experi ments 
should constrain precisely (see e.g. 
Eiscnstcin et al.| 1999). 



Jungman et al. 1996; 



As an example, we take a sample with z mo dian = 1 and 
n — 2 x 10 5 de g~ 2 as appro priate for a magnitude limit 
of R ~ 25 (see Smail et al.| 1995b). The signal to noise 
in the full sample is quite high, e.g. at t = 100, S/N ~ 
25. Thus we expect that subdividing the sample should 
improve parameter estimation. 

TABLE 1 

Parameter Estimation for z mcd ian = 1 

1 ~When taking these derivatives the redshift distribution 7lj(z) is 
held fixed as opposed to the distance distribution [ni(z)dz/dD] in 
equation (ph. 
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As shown in Table 1, subdividing this sample in equal 
halves, denoted as 2(1/2), improves the errors a a by a 
factor of 2 to 7. Since the signal in the lower redshift 
bin is smaller than in the higher redshift bin, it suffers 
comparatively more from the intrinsic noise variance. One 
can optimize the binning to correct for this effect. Divid- 
ing the sample so as to isolate the upper quarter [2(1/4)] 
improves the errors modestly whereas isolating the upper 
eighth deproves them. We plot the full range as a function 
of the fraction of galaxies in the upper bin in Fig. |^. Notice 
that though the improvement factor is roughly flat from 
0.15 — 0.5, it drops rapidly when noise dominates either 
the upper or lower fraction. If the signal were the same in 
both bins, this would occur at 0.04 and 0.96 for I — 100. 
The fact that the true improvement is skewed to smaller 
upper fractions reflects the fact that the signal increases 
to higher redshifts. 

Moving to three divisions makes only a small improve- 
ment over two. In Table 1 we give the results of taking 
3 bins with an equal number of galaxies in the upper two 
bins, e.g. [3(1/4)] represents a division by number of (1/2, 
1/4, 1/4). In fact the errors for three bins can be higher 
than those with two if not chosen wisely. 

We conclude that for a redshift distribution of the form 
given by equation (|^) with 2 me dian = 1, a = 1 and /3 = 4, 
crude partitioning suffices to regain most of the redshift 
information in adiabatic CDM models where the change 
in the growth rate across the distribution is slow and con- 
trolled by a small number of cosmological parameters 

How robust are these conclusions against changes in the 
distribution and model? A wider redshift distribution of- 
fers greater opportunities for tomography. For example, 
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let us widen the distribution by taking (3 = 2 in equa- 
tion (§). Then the gains by simply halving the distribu- 
tion are a factor of 9.7 for J7a; going to a 3(1/4) scheme 
raises this to 12. 

These considerations are also relevant for deeper sur- 
veys. Consider a survey with z mc dian = 2 and n = 3.6 x 10 5 
deg -2 . The parameter estimation results are given in Ta- 
ble 2. Not only is the overall improvement from subdivi- 
sion larger (up to a factor of 24 for three bins) but the 
relative improvements between parameters also changes. 
This is because even within the adiabatic CDM paradigm 
the importance of the different parameters in determining 
the growth of structure depends on rcdshift. 

Perhaps the most important aspect of weak lensing to- 
mography is that it has the ability to falsify the underlying 
adiabatic CDM model. For this reason, it is wise to ex- 



amine the power spectra from the rcdshift bins directly, 
since these are the observables, rather than jump directly 
to modelling the data with parameters under the adia- 
batic CDM framework. For example, tomography may 
show that the component that accelerates the expansion 
of the universe is not the cosmological constant or call into 
question the fundamental assumption that structure forms 
through the gravitational instability of cold dark matter. 

TABLE 2 

Parameter Estimation for z mo dian = 2 
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4. DISCUSSION 

We have shown the precision with which cosmological 
parameters can be measured from a weak-lensing survey 
can be significantly enhanced by tomographically deter- 
mining the evolution of the statistical properties of large- 
scale structure across the finite redshift width of the source 
distribution. Crude redshift binning of the data can re- 
cover most of the statistical information contained in the 
redshifts. For example, most of the gain for a magnitude 
limited survey with z me( j; an = 1, under the adiabatic cold 
dark matter paradigm, comes from separating out the up- 
per and lower redshift halves of the distribution. For wider 
distributions and stronger rates of change in the growth of 
structure, more information can be extracted by finer bin- 
ning, especially of the higher redshift portion of the sam- 
ple where the signal is greater. The appropriate number 
of bins can be empirically determined by examining the 
correlation between bins and the noise properties of the 
data. 

We have been assuming that the individual redshifts of 
the galaxies will be known sufficiently precisely to deter- 
mine the redshift distribution of the subsamples. Realis- 
tically, the redshift information will be limited by the ac- 
curacy of photometric redshift techniques which cur rently 
show errors of Az - 0.1 (68% CL) for 0.4 £ z £ 1.4 ( [Hogg 
1998). While statistical errors on the large samples 



et al. 



rors or biases in the technique may cause problems. It is 
beyond the scope of this letter to test these issues fully. 
To give some feel for their effect, let us consider the me- 
dian redshift z mc dian as an additional parameter with a 
prior uncertainty from photometric redshifts of the full in- 
dividtual error 0.1. Including this uncertainty degrades 
the precision in the parameters by 3% in the worst case. 

While this effect is negligible, more worrying is a bias 
that is a function of redshift, especially in the largely 
untested regime 1.4 <J z <; 2, as that can shift the dif- 
ference between the power spectra of the subdivisions. 
Isolating the few percent of galaxies at z > 2.5, where 
the techniques are tested, yields gains that are compara- 
ble to the optimal division (see Fig. || lower panel), but 
the compactness of such galaxies poses an obstacl e for 



measuring the lensing distortion from the ground (Stci- 
del et al.| 1996) 



hows 



Despite these caveats, this study s : 
that tomography with weak lensing is both possible and 
would substantially improve the precision with which we 
can measure the growth of structure in the universe. 
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